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parton distributions for the p — > At (and p — > A £ ) transition and a hard subprocess amplitude for 
ww — > cc. Thereby the mass of the charm quark is taken as the hard scale so that our results are not 
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parton distributions for the p — > Ac transition by an overlap of simple quark-diquark light-cone 
wave functions we make first predictions for p^A c transition form factors and unpolarized pp — > 
A £ Ac cross sections. Our findings may become interesting in view of forthcoming experiments at 
FAIR in Darmstadt. 
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1. Introduction 

Exclusive hadronic reactions which require the production of heavy quark-antiquark pairs are, 
to a large extent, cleaner and easier to handle than those in which only light flavors are involved. On 
the one hand, certain elementary reaction mechanisms can already be ruled out from the beginning 
due to the small heavy-flavor content of the quark sea. On the other hand, the mass of the heavy 
quark itself can serve as a hard scale so that there is a good chance that QCD perturbation theory 
provides a substantial part of the process amplitude already at moderately large energies. 

The simplest elementary reaction mechanism for the process we are interested in, namely 
pp — > AfAf, is depicted in Fig. This is the mechanism which we assume to be dominant in the 
forward hemisphere for energies well above production threshold. We are going to analyze it in 
terms of generalized parton distributions [[]]], [||]. 




Figure 1: Double-handbag contribution to pp — > A C A 6 . 



2. Kinematics 

We denote particle momenta and helicities as shown in Fig. |l|. For our purposes it is most 
convenient to work in a center-of-mass frame, in which the light-cone (LC) components of the 
proton and A c momenta are parameterized as 1 



m 2 + A{/4 
2(1 + £)P + 



' 2 



and 



M 2 + A 2 L /4 Aj_ 
2(l-£)p+'~T 



(2.1) 



respectively. M stands for the A c mass and m for the mass of the proton. The corresponding 
antiparticle momenta are q = [p~ ,p + , Aj_/2] and q' = \p'~,p' + , — Aj_/2]. The average of proton 
and Ac momenta p = \ {p + p') defines the longitudinal direction and the 4-momentum transfer is 
specified by A = p' — p = q — q'. The relative momentum transfer in longitudinal direction is given 
by the "skewness parameter" 

- " + -^ 4+ (2.2) 
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We use the following convention for the LC components of a 4-vector v: 



,v ,Vj_ 



withv* = 1/V2(v° ±v 3 



and v 
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3. Factorization 

If the dynamical mechanism underlying pp — * A C A C scattering is given by Fig. [I], the corre- 
sponding scattering amplitude can be written as {kf = + &-)/2) 



M , i 
fi v ,nv 



Jd 4 k?e(kT + ) f±*&«j*l%e{$-) J^e ik ^H(k[,k' 2 ;k u k 2 ) 



x (A, : p 4 i \TW ( -|) ¥" (|) \p : p^)(A c : q ',v'\TV" (|) «F e (-|) |/> : </. v) 



(3.1) 



with H (k\,k' 2 ;ki,k2) representing the scattering amplitude for the elementary subprocess 

u(k\, Ai) u{k2,^i) - ► Aj)c(^ 2 , A2) . (3.2) 



The hadronic matrix elements describe the emission of a light (anti)quark by the (anti)proton and 
the absorption of a charm (anti)quark by the (anti)A c . For better readability we have suppressed 
helicity labels as well as color and spinor indices (and corresponding sums) for the quarks. 



In order to simplify the right-hand-side of Eq. ( p.lD we make a few plausible assumptions. 
In analogy to Compton scattering [||], [Q] we assume that the process is dominated by partons 
with restricted virtualities and transverse momenta. To be more precise, we demand for the ac- 
tive (anti)quarks that kf < A 2 and \k\ — nr c \ < A (cf. Fig. [l] for the assignement of (anti)quark 
momenta). The relative transverse momentum components should satisfy 2 

k 2 ±i / Xi <A 2 , £?,./*/< A 2 , (3.3) 

with A being a hadronic scale of the order of 1 GeV and m c the charm-quark mass. The restrictions 
for the transverse momenta are supposed to hold in the hadron in/out frames (tilde/hat), in which 
the corresponding parent hadrons carry no transverse momenta. Analogous constraints should also 
hold for spectator partons. Furthermore, the hadronic matrix elements represented by the blobs 
in Fig. [l] are assumed to exhibit a pronounced peak at xq = m c /M at 0.6 — 0.7. This assumption 
reflects the expected behavior of light-cone wave functions for the A c and is also supported by 
similar observations for heavy-quark fragmentation functions. 

Since the virtuality of the gluon in the subprocess amplitude has to be at least 4m 2 , it occurs 
to be natural to take the charm-quark mass m c as a hard scale. Under the foregoing assumptions 
on the parton momenta it is then a good approximation to neglect the relative transverse momenta 
(kxi and k' ±i ) of the (anti)quarks in the subprocess amplitude H (k[,k' 2 ;ki,k2) and to replace the 
momenta of the active quarks kf by (on-shell) momenta kf^ which are collinear to those of the 
corresponding hadrons. With these approximations in the subprocess amplitude most of the inte- 
grations in Eq. ( |3~l| ) can be done analytically. The resulting delta functions enforce a light-like 
separation of the fields in the hadronic matrix elements. This allows us to drop the time ordering 
of the fields. Due to our assumption on the hadronic matrix elements it is also justified to apply a 



The plus-components of parton momentum k and hadron momentum p are related by k + = xp + . For further 
purposes it is also convenient to introduce an "average" momentum fraction of the active quarks, i.e. x\ = (£[ V+ )/ (P + ) 
(and analogously for the active antiquarks) 
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"peaking approximation" to the subprocess amplitude (i.e. replace all x\ by xo). In this way we are 
able to write the hadronic amplitude as a product of integrals over soft hadronic matrix elements 
with a hard partonic subprocess amplitude: 



x [dk™ + d(k\ v+ 



2% 



\P ■ P,W 



J- 



dz 



x / d*f-0(*f ") / ^-e ik > * (Ac : q,v\¥ u I ) V c I ) |p : g> V 



27T 



^2 



(3.4) 



The subprocess amplitude can be calculated perturbatively, whereas the hadronic matrix elements 
describe the non-perturbative transitions of the (anti-)proton to the (anti-)A c by emission of a (anti)w 
quark and reabsorption of a (anti)c quark. It are just the hadronic matrix elements which give rise 
to generalized parton distribution functions (GPDs) and to form factors, which are essentially 1 jx 
moments of the GPDs |JIjl, Therefore we will have a closer look on them. 

After some algebraic manipulations (insertion of energy- and helicity-projectors, etc.) the first 



double integral in Eq. ( p.4[ ) can be split up into three (independent) terms of the form 

d ^= p + [ ^Le*^ (Ac : pit' |¥ c r¥" \p : pa) =: f ^ 
2n ' 



1 



jr?(r), (3.5) 



fx n 



where T = y + , y + Js, or ic +] . y + and Y + Y$ demand for the same helicity of the emitted and absorbed 
parton, whereas /a +J is associated with a helicity flip. At this point the GPDs are introduced by 
decomposing M"f (r) into its covariant structures: 



^(r + r 5 ) = u[p',u' 



^■ r ;;(i° +j ) = u(p,u' 



H™(x,t;,t)r + +E c «(x,t;,t] 



/C7 +V A V 



H cu (x,^t)r + Y5+E cu (x^,t 



M + m 

y 5 A+ 



u(p,jl) 



M + m 



u(p,u) 



H? (x, $,t) ia + J + H T U (x, t;,t) P -^ z P -+ 



E T " (x, § ,r) Y+A ' ■ A+f + E C T U (x, Z,t) 



Mm 

Y^pj — p + Y j 



(3.6) 
(3.7) 

(3.8) 



u(p,u). 



M + m ' 1 (M + m)/2 

Analogous results hold for the hadronic matrix elements of the antiparticles. 

What enters the scattering amplitude are rather moments of the GPDs than the GPDs them- 
selves. It is thus convenient to introduce form factors 



dx\ 



(3.9) 



where F, stands for any of the transition GPDs introduced in Eqs. (3.6) - (3.8). In our notation the 
form factors Ry, Rt, Ra, Rp, St, Syi, Ss, Sy2 correspond to the GPDs H, E, H, E, Hj, Et, Ht, Et, 
respectively. 
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Combining Eqs. (3.5) - (3.9) and inserting them into Eq. (3.4) allows us to write the hadronic 
scattering amplitude in terms of these, a priori unknown, form factors and the hard subprocess 
amplitude. The result is a rather lengthy expression. Fortunately, some of the form factors can be 



neglected. By taking different helicity combinations for the mauix elements in Eqs. ( p.q ) - ( |3.8| ) 
we observe that all the form factors apart of Rv,Ra and St (corresponding to the GPDs H,H and 
Hj) must vanish in forward direction, since they are connected with orbital angular momentum. 
Also for large Mandelstam t these form factors should be suppressed (as compared to Rv,Ra and 
St) by powers of \f—i, depending on how many units of angular momentum are involved. We 
therefore expect the following hierarchy of the form factors: 

\Rv\,\R A \,\S T \ » \R t \,\Rp\,\S s \,\S vi \,\S V2 \. (3.10) 

In M^iyi ^ v (p' ,q'; p, q) a further suppression of contributions coming from helicity flip form factors 
is caused by the corresponding subprocess amplitudes. We therefore neglect all helicity-flip form 
factors apart of St- This provides us with comparably simple expressions for the pp — > A C A C 
helicity amplitudes in terms of the form factors Rv,Ra and St and the subprocess amplitudes: 

M+± ' +± = ^C 1- £ 2 ) (4 =F*a), (3 ' U) 
M—, ++ = 0, (3.12) 
M-+,+- = ^(l-£ 2 )//- + .+-4, (3.13) 

M ++ ,± T = ^-{1-^ 2 )H ++ . + ^St(R v +Ra), (3.14) 



c 



M T ±, ++ = t^{1~?; 2 )H ++ , + -St(Rv-Ra), (3.15) 



with 



#+_,+_ = -Siias(4s) cos 2 (0 cm /2) , #_+,+_ = Sneers) sin 2 (0 cm /2) , (3.16) 
#++,+_ = -Sn:a s (4s)-r sinQcm, H++ t -+ = H+++- . (3.17) 

Cp =4/3 denotes the color factor, Nc = 3 is the number of colors. The remaining helicity ampli- 
tudes are related by parity invariance, 

M-p-y-n-v = (-l) M '~ v '^ +v M AVv . (3.18) 

For a more detailed account of these results and their derivation we refer to Ref. [§]. 

4. Modelling the GPDs 

The kinematical requirement for the production of a cc pair x > 2xqM /yfs implies that x > E, ? 
We are thus in the DGLAP region where GPDs can be modelled as overlap of light-cone wave 
functions of the proton and the A c [Bj. For the A c it is certainly a good approximation to consider 



This is also the reason why the integral in Eq. (3.9) starts from t, and not from 0. 
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only its valence Fock state. For the proton higher Fock states may be important, but the required 
overlap with the A c projects out only appropriate spin-flavor combinations of its valence Fock state. 

As a first attempt we consider a simple quark-diquark model for the baryons. We assume that 
a baryon consists of an active quark, which undergoes the hard scattering, and a diquark, which 
acts as a spectator. For the A c the diquark has to be a spin-isospin scalar. In Ref. [^] pp — ► A C A C 
has already been studied within a quark-diquark model, but without using the general framework of 
GPDs. For our present investigation we will take a quark-diquark wave function for the A c which 
is similar to the one used for a calculation of heavy-baryon transition form factors [|l] : 



|A,1 



\c±S\ 



nil] , 



with x ¥ A (x,k x )=NA{\ 



-x)e~^=3 K ^e~ a * M -&=V. (4.1) 



x denotes the momentum fraction carried by the quark. The quark-diquark wave function of the 
proton, on the other hand, is chosen similar to the one in Ref. [Q] : 

a- 2 

\p:±)(x,k x ) = \u±S [udl ) with W p (x,kx)=N p (l-x)e-*&l kx . (4.2) 



These are pure s-wave wave functions and hence all GPDs apart of H, H and Hj will vanish. 
Within this model the quark should have the same helicity as the baryon, which implies H = H. 
Likewise, also Hj is non-zero and we even have 



Ht — H — H , 



(4.3) 



i.e. we are left with a single GPD. For the model wave functions ( J4- 1| > and ( |4.2D the overlap integral 
can be done analytically with the result 



1 



;i-x) 3 (*- 2 -^) 



x exp 



(N A N P \ 



-a K M 



exp 



a\a\{\ 



al(l-^ 2 (x + ^+aj(l + ^ 2 (x-^ 



(4.4) 



With this GPD we can now evaluate the form factor R = Rv=Ra = St- 



5. Numerical Results and Conclusions 

For our numerical calculations the transverse size parameters a A = a p = 1 . 1 GeV~ 1 are chosen 
in such a way that they provide a value of « 300 MeV for the mean intrinsic transverse momentum 
< k\ > x l 2 of a quark inside a baryon. The wave function normalizations are fixed such that the 
probabilities to find the quark-diquark states in a proton or a A c are 0.5 and 0.9, respectively. The 
x-dependence for the resulting transition GPD H at different values of t (i.e. different values of 
the skewness parameter is shown in Fig. ^[ Shown in the same figure is also the corresponding 
(scaled) transition form factor 7? as a function of t. With this form factor the unpolarized differential 
and integrated cross sections for — > A f A c are easily evaluated. Corresponding predictions are 
plotted in Fig. SI. The integrated cross section is of the order of nb, i.e. comparable in magnitude 
with the result in Ref. [^]. This is still in the range of high precision experiments. Whether pre- 
dictions from more realistic three-quark wave functions will be comparable to the outcome of this 
simple quark-diquark model is presently under investigation 
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Figure 2: Le/f: The transition GPD H in the quark-diquark model versus x at Mandelstam s = 
30 GeV 2 for A]_= (solid), 2 (dashed), 4 (dotted) GeV 2 (corresponding to -t = 1.13, 3.30, 5.54GeV 2 or 
E, =0.11,0.12,0.14). Right: The transition form factor R scaled by \t\ in the quark-diquark model versus 
\t\ [GeV 2 ] for Mandelstam s = 30 GeV 2 . 




Figure 3: Left: The pp — > A C A C differential cross section do/dQ. versus cos0 C m at Mandelstam s = 
30 GeV 2 . Right: The pp ^ A C A C integrated cross section c versus Mandelstam s. 
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